Schwarzschild black holes can wear scalar wigs 
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We study the evolution of a massive scalar field surrounding a Schwarzschild black hole and find 
configurations that can survive for arbitrarily long times, provided the black hole or the scalar field 
mass is small enough. In particular, both ultra-light scalar field dark matter around supermassive 
black holes and axion-like scalar fields around primordial black holes can survive for cosmological 
times. Moreover, these results are quite generic, in the sense that fairly arbitrary initial data evolves, 
at late times, as a combination of those long-lived configurations. 
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I- Introduction. It has been known for some time that a 
Schwarzschild black hole (BH) cannot support a nontriv- 
ial scalar field (SF) distribution, i.e. such a hole does not 
admit scalar hair [HQ. However, the no-hair theorems 
do not exclude the existence of dynamical solutions that 
decay very slowly in time. In a recent work Q we found 
regular SF configurations surrounding a Schwarzschild 
BH that can survive for relatively long times, and con- 
jectured that, for a certain range of values of the SF and 
BH masses, such SF configurations could survive in the 
vicinity of the BH for cosmological timescales. Here we 
show that this is indeed what happens, and that such 
distributions arise from the evolution of fairly arbitrary 
initial configurations. Our results give support to the 
idea that the dark matter halos could be described by a 
coherent scalar excitation 0,(1 (see for papers with 
a similar motivation) , and to the possibility of nontrivial 
axion distributions surrounding primordial BHs 0. 

II- Three types of resonances. The dynamics of a 
SF propagating on a Schwarzschild background are de- 
scribed by the Klein- Gordon equation 



d 2 <t> 



dr 2 



V{r)4> = 



with V(r) an effective potential given by 
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Here M is the BH mass, £ the angular momentum number 
(we decompose the field in spherical harmonics), yU the 
mass of the SF (in units for which G = c = H = 1), and 
r* (— oo < r* < +oo) the radial "tortoise" coordinate, 
related to the Schwarzschild radial coordinate r (r > 2M) 
by = r + 2Mln(r/2M — 1). As discussed in detail 



in Q (see also reference @), for each £ > and Mfj, < 
M/j, CI it(£), V describes a potential well which is enclosed 
between a barrier close to the BH, and the asymptotic 
positive value [i 2 of the potential at infinity. 

The stationary solutions of Eq. (JlJ are characterized 
by a time dependency of the form e %u)t with a real fre- 
quency uj. Such solutions behave as a combination of 
ingoing and outgoing waves to the left of the potential 
barrier, close to the horizon, and we require them to de- 
cay exponentially to zero at spatial infinity. There is a 
continuous spectrum of such solutions for to 2 < fj, . Fur- 
thermore, as we have shown in for to lying in the 
region V nlln < oj 2 < /i 2 , with V m - ln the local minimum of 
the potential, there is a discrete set of such frequencies for 
which the amplitude inside the potential well takes very 
large values when compared with the amplitude close to 
the horizon. We call these modes the stationary reso- 
nances, and the region V m - m < oj 2 < fj? the resonance 
band. 



If one imposes the condition of no waves coming from 
the region close to the horizon, while keeping the require- 
ment of exponential decay at spatial infinity, it turns out 
that Eq. (TTJ) is only satisfied for a discrete set of com- 
plex frequencies. These solutions have been called quasi- 
resonances in the literature [Io| (see also 11 1 for a pre- 
vious study of such solutions). 

Both the stationary resonances and quasi-resonances 
are in fact non-physical solutions. First, as mentioned 
above, all purely stationary solutions, i.e. those with 
real u, require waves to move outward from the hori- 
zon region to compensate for the waves that tunnel out 
through the barrier and move toward the horizon (as oth- 
erwise the situation would not be stationary). Imposing 
the condition of no waves coming out from the horizon 
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clearly improves the situation at the cost of introducing 
complex valued frequencies u. This makes sense physi- 
cally since now the solutions must decay in time as the 
waves tunnel out of the potential well and fall toward the 
BH, with this decay represented by the imaginary part 
of the frequency. Nevertheless, it turns out that such so- 
lutions are still non-physical, as one can show that the 
energy density diverges at the horizon for both types of 
solutions. 

There is yet a third class of resonant solutions we will 
consider. In reference 0] we performed numerical evolu- 
tions of regular initial data that has not such divergence 
of the energy density at the horizon, and found damped 
oscillating solutions that remain surrounding the BH for 
long times. These solutions have frequency equal (within 
the numerical error) to that of the stationary resonances. 
We call these solutions dynamical resonances. We will 
show that the real part of the frequency of the quasi- 
resonant modes coincides with the frequency of oscilla- 
tion of the stationary and dynamical resonances, and that 
the imaginary part coincides with the decay rate of the 
dynamical resonances. This is a nontrivial result since, 
after all, the dynamical solutions arc in some sense "in- 
finitely different" to both the stationary resonances and 
the quasi-resonances because they are regular solutions 
with finite energy. This is similar to the case of the quasi- 
normal modes of BHs, which also diverge at the horizon, 
but nevertheless physical excitations behave as combina- 
tions of them locally at late times. 

777- Toy model. In order to illustrate how the different 
resonant modes arise and are related to each other, we 
first consider a simple toy model in which the potential 
V in Eq. (|lb[) is replaced by 

U(x) = A5(x) + p?0(x- a) , x := r* , (2) 

with A and a positive constants. That is, we have a 
potential well enclosed between the <5-barrier and a step 
function where U jumps from zero to its asymptotic value 
p 2 . For Mp <C 1 and I > 1 the minimum of the 
potential well of the physical potential V is located at 
r m in = £(e+l){AIp 2 )- 1 [l + 0(Mp) 2 }, and the strength 
of its barrier is j^™ n Vdr* = [U{1 + 1) + l](4A/)- 1 [l + 
O(Mp) 2 }. Therefore, we require the dimensionless pa- 
rameters A := 2p/A and p := (a/i) _1 to scale like Mp 
in our toy model potential U. In the analysis below we 
prove that for fixed p < 1/n and small enough values 
of A, the potential U gives rise to a discrete set of sta- 
tionary and quasi-resonant frequencies whose real parts 
agree with each other. 

The resonant modes for the potential U have the form 



(f)(t,x) 



ae luJX + (3e~ lu ' x , x<0, 
e-vV-" 2 * x > a , 



(3) 



with uj a complex frequency. We require u (— oo, — p] U 



therefore we can choose the sign of the square root such 
that its real part is positive, leading to exponential de- 
cay for x > a, as required. The constants a, (3, 7 and 
S are determined by the matching conditions at x = 
and x = a which consist in the continuity of 4> and its 
first spatial derivative <p' at x — a and the continuity of 
</) and the jump condition tf>'(t, 0+) - <j>'(t, 0~) = A<p(t, 0) 
at x = 0. 

For the stationary resonances to is real and our restric- 
tions imply cj G (— p, p). Since U is real, it follows in this 
case that a = /3*, 7 = 5*, and the ratio between the am- 
plitudes of the solution for x < and the one inside the 
well is v(u>) = |a|/|7|. The resonance frequencies are the 
ones that minimize this ratio. For the quasi-resonances, 
in turn, the condition of no waves coming from the hori- 
zon implies that j3 = 0, which can only occur for a non- 
trivial imaginary part of u. In terms of the quantities A 
and p defined below Eq. ((2|), and the complex angle tp 
lying inside the strip — n/2 < Re(ip) < n/2 defined by 
lu = /isiniy2, the condition /3 = reads 



_F p (A, tp) := 1 + iA sin ip — e 



—2itp—— sin ip 



0. 



(4) 



while the ratio of interior and exterior solutions in the 
stationary case is 

1 \F P {\M\ 
A sin ip 



(5) 



As we prove now there exists, for each p < 1/n and small 
enough A, a finite family <p n (A) of solutions of Eq. Q 
with the following properties: they depend smoothly on 
A, their imaginary parts decay like A 2 , and their real 
parts agree with the local minima of v(u>). This shows 
the existence of quasi-resonant modes for the toy model 
potential U whose frequencies are related to those of the 
stationary resonant modes through their real parts. 

In order to show this we consider first the limit 
A = 0, in which case Eq. yields real solutions tpn G 
(— 7r/2,7r/2) satisfying simp^ = p(nix — ipn ), with 
n = 1,2,3,.... There are a finite number of solutions 
when p < l/?r. For p -C 1/V the fundamental solution 
is approximately tp^ ~ p-K. In order to extend these 
solutions for small A > we invoke the implicit function 
theorem and observe that the function F p (X,p) is lin- 
ear in A and analytic in p. Furthermore, (0, ) = 

2i{l+p- 1 cosp i n ) ) ^ since -n/2 < p { n ] < n/2. There- 
fore, there exists for small enough A a unique smooth 
solution curve tp n (\) satisfying F p (X, ip n (X)) = and 

fn (0) = <Pn^ ■ A Taylor expansion reveals that the corre- 
sponding complex frequencies w n (A) = psn\{p n {\)) have 
the form 



Re(cj„(A)) = ps\np ( n 



(0) 



1-- + ^X 2 +0(X 3 ) 
In qi 



[p,, 00), which guarantees that Rey/T 2 



7^ 0, and 



• 2 (0) 

ImK(A)) = ^^L^A 2 +0(A 3 ) 



,(6a) 
(6b) 
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where we have defined q n := 2 (cos 1 ip^ + p 1 ) and 
a := 1 — q^ 1 sin 2 <pn / cos 3 ifn . Furthermore, one can 
verify that the local minima of the function v{uj) given 
in Eq. ([5]), describing the stationary resonances, have pre- 
cisely the expansion given in Eq. (|6a[) . Therefore, up to 
second order in A, the frequencies of the stationary res- 
onances agree with the real part of the quasi-resonant 
frequencies. 

Since both parameters A and p scale like Mfi in our 
model potential, it follows that for Mfi <C 1 the imagi- 
nary part of Mu> n scales like (Mfi) 6 . This explains, at 
least qualitatively, why the quasi-resonant modes are as- 
sociated with such large time-scales when Mfi <C 1. 

It is also possible to relate the quasi-resonant frequen- 
cies uj n to the dynamical resonances by representing the 
solutions of Eq. ([T]) as an inverse Laplace integral. This 
calculation will be presented elsewhere. 

IV- Frequency comparison and characteristic times. 
For the physical potential, Eq. (jlbjl. the quasi- resonant 
modes can be obtained semi-analytically in several ways, 
the most common are the continued fraction method in- 
troduced by Leaver [l2[ , and the WKB approach (see 
l44l7| for details) . An analytic expression valid for small 
values of the combination Mfi was obtained by Detweiler 
in [l0 | . The stationary and dynamical modes were previ- 
ously considered in [3[ . In this section we show that the 
oscillations of the dynamical resonances have frequencies 
equal (up to numerical error) to the real part of the fre- 
quencies of the quasi-resonant modes, while the decay 
rates are equal to the imaginary part of such frequen- 
cies (see Fig. []}. We also show that the imaginary part 
of the quasi-resonant frequencies calculated by means of 
Leaver's method has the limit behaviour predicted by 
Detweiler. Hence, we can use Detweiler expressions to 
determine the decay rate of the dynamical resonances 
with very small values of M ft,, which cannot be reached 
with the other two methods. 

Leaver's semi-analytical approach [l2[ assumes a har- 
monic time dependence for the SF, 4>(t,r) = ij)(r)e lu)t , 
with uj the complex quasi-resonant frequency, and pro- 
poses a power series expansion for the radial function 
ip(r) of the form 



, ^ / 2M\ n 

tp = V > J On [ 1 - — I , 



n=0 



(7) 



where, following 



Ref. [181, 

2M/r) 

antees that close to the horizon ip 



we 

2iMu 



define ip := 
This guar- 
e %UJT , whereas at 



spatial infinity -0 



, with x = Vm 2 



and the sign convention for the square root was ex- 
plained in the previous section. Notice however that our 
boundary conditions are different than those in [l8j |. 

In Fig. [T] we plot the imaginary part of the frequency 
of the first quasi-resonant mode (the one with the lower 



value of Re(Mw)) for the cases £ — 0, 1 and different val- 
ues of Mfi. We find that the real part of the frequencies 
obtained corresponds to the frequencies at which the SF 
oscillates in a dynamical scenario Q. Furthermore, we 
also find that the imaginary part of the frequencies (dots 
in Fig.Q]) corresponds to the decay rate of the field in the 
numerical evolutions (empty circles in Fig. [I}. 

In the numerical evolutions we were not able to evolve 
the configurations for small values of Mfi, both because 
of the propagation of the errors in the numerical approxi- 
mation, and because the time required for such evolutions 
becomes prohibitive. Even though Leaver's method al- 
lows us to obtain results for parameter values that could 
not be reached with the numerical evolutions, numerical 
roundoff errors make it still prohibitive to obtain accu- 
rately very small values for the imaginary part of the 
quasi-resonant frequencies. 
In reference 
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analytic expressions for the spectrum 
of quasi-resonant modes were found in the limit Mfi <C 1, 
for values of £ > 1. For I = 1, the imaginary part of the 
frequency for the first quasi-resonant mode in this limit 
is given by Im(Mw) = (Mfi) w /6. Strictly speaking, the 
expressions given in (llj may not apply for the case with 
£ = 0. However, if we nevertheless use them in that case 
we find Im(Muj) = 16(M/i) 6 . We have calculated the 
quasi-resonant mode frequencies using Leaver's method 
up to the smallest values of Mfi allowed by the round- 
off errors, and we have found that the imaginary part of 
such frequencies matches the Detweiler approximations 
(within errors) at small values of Mfi (see Fig. [T]). This 
result allows us to conclude that, for small values of Mfi, 
the decay time of the dynamical resonances should cor- 
respond to the one predicted by Detweiler in (llj . 




FIG. 1: The imaginary part of the frequency for the first 
quasi-resonant mode obtained using Leaver's method is shown 
as function of Mfi for the cases I — 0, 1. 

Since the decay rate of the dynamical resonances is 
related to the imaginary part of the quasi-resonant fre- 
quencies, then its half-life time is inversely proportional 
to Im(Muj). In particular we can find BH and SF masses 
such that ti/2 ~ 10 10 yrs. and larger for the cases £ = 0, 1 
previously mentioned. Those results are shown in Fig. [21 
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There are two distinct regions of the parameter space of 
physical interest for which the configurations live longer 
than the age of the Universe: (i) a SF mass smaller than 
1 eV and BH mass smaller than 10~ 17 M Qj consistent 
with primordial BHs with an axion distribution Q; and 
(ii) an ultra-light (fuzzy) SF [Hil with mass smaller 
than 10 -22 eV and a supermassive BH with mass smaller 
than 5 x 1O 1O M , as could be the case for a dark matter 
halo [1, [H surrounding a BH at a galactic center. 




t„ 2 =10 yrs(l=I) - 
t 1/2 =10 10 yrs(l=0) 
t |;1 >10 yrs for all 1 



FIG. 2: Combinations of the parameters M and fi for which 
the SF half- life time can be larger than the age of the universe. 

V- Evolution of arbitrary configurations. In this 
section we study the long-term numerical evolution 
of SF distributions that are initially surrounding a 
Schwarzschild BH. We show that at late times, after 
some SF falls into the BH and some escapes to infin- 
ity, the remaining SF consists of a superposition of the 
dynamical resonances, which, for practical purposes, can 
be described in terms of the quasi-resonant modes. Re- 
markably, we see that this holds even for pretty arbitrary 
initial data. Note that the calculations in this section are 
identical to those in but with the crucial difference 
that here they are performed on arbitrary initial data 
(as opposed to initial data obtained from stationary res- 
onances), hence leading to new, much stronger results. 

In order to test the evolution of a variety of SF distri- 
butions, we construct a two-parameter family of initial 
data, of the form 



Uo(r) = 



N(r - i?i) 4 (r - R 2 y for R 1 < r < R 2 
I otherwise 

(8) 

with the normalization TV = [2/(i? x — i? 2 )] 8 . The parame- 
ters i?i and i? 2 arc chosen in order to set distributions at 
t = with different "locations" and "sizes". The initial 
value of the time derivative, u\ t=Q , is constructed as 

, f(t,r) 



dj_ 

dt 



uo(r — vt) 



(9) 



figurations studied is simply set to zero. The evolution 
equation is solved numerically, as described in [22| . 

As expected, during the initial stages of the evolution 
some SF accretes into the BH, while some is radiated 
away, both with rates that depend greatly on the initial 
data chosen, and can be very large in some cases. How- 
ever, at late times, all evolutions show a similar steady 
behaviour with slow accretion into the BH. Similar re- 
sults arc obtained when studying the long-term evolution 
of a variety of configurations, some of them very differ- 
ent in size and spatial distribution —although all distri- 
butions studied here are quite "wide" when compared to 
the BH size. 

One can gain some understanding of this behaviour by 
conducting a spectral analysis as follows. We perform a 
(discrete) Fourier transform in time, ^[(^(t, rj)](oj), of the 
SF at several (about 20) sample points rj, then take the 
average to obtain the spectrum J r [(/)(t)](w). For brevity, 
we present here results corresponding to only two rep- 
resentative examples. They correspond to initial data 
with i?i = 100A/, i? 2 = 200AZ; and i?i = -I000M, 
i? 2 = 1000A/. In both cases (mM = 0.3, I = 1 and v = 0. 
The spectra are shown in Fig. [3l 




where v is a free parameter that has no significant effect 
on the results presented here, and for the majority of con- 
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FIG. 3: Fourier transform in time of the evolution of initial 
data with i?i = 100, R 2 = 200 (top panel); and i?i = -1000, 
i?2 = 1000 (bottom panel). The peaks to the right of the res- 
onance band in the top panel are a numerical artifact caused 
by noise originating at the outer boundary. 

For comparison, we have indicated in the figure the 
first resonance frequencies (vertical lines), as obtained 
in One can clearly see peaks that coincide with each 
resonant frequency. These results seem to indicate that, 
at late times, even quite generic SF distributions evolve 
as a combination of the resonant modes, which, as we 
have shown, can last for cosmological time-scales. 
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